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Most of the common used models of epidemic spreading allow contaminating many neighbors of 
a particular node in the network. They are usually analyzed by differential equations on probability 
vectors. We propose a model of epidemic spreading, which restricts to at most one contamination 
per time step and analyze it by discrete approach, working on vectors of possible states of the 
system. Theoretical predictions of epidemic treshold, stationary state and time needed to reach it 
are given and appear to be perfectly consistent with computer simulations. We also point ou that 
the model appears to be well suited to mimic epidemic spreading within student communes. 
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I. INTRODUCTION 

It has been admitted, that the most appropiate m.od- 
els of epidemics spreading are these based on dynamical 
processes on particular graph models rather than those 
defined by fenomenological differential equations [l|, y]- 
Within this approach the nodes of a network are usu- 
ally considered as individuals, who are connected with 
each other by vertices corresponding to social links. Al- 
though some authors use continous time simulations (see 
G-g- 3)i the approach presented commonly (see [J| for 
a review) is based on the idea that at each discrete time 
step a particular node of the network can contaminate 
each of its neighbors with some finite probability p. The 
whole set of vertices is being divided into compartments, 
usually referred to susceptible (S), infected (I) and recov- 
ered (R) individuals, but the general mechanism stays 
more or less unchanged. There has been a broad range 
of methods developed in order to analyze such models. 
In the most basic approach people assume individuals 
to be identical and homogeneously mixed (homogeneous 
assumption, [4|). In order to take into account hetero- 
genity of the system a kind of block approximation hax 
been used Q, treating nodes with the same degree as 
statistically equivalent. This is not always enough, as 
some real networks manifest degree correlation, mainly: 
the conditional probability, that two vertices of degree 
k, k' are connected depends on both degrees k, k' p. 
The next step thus is to take into account correlation 
[7|. Finally, one can employ whole adjacency matrix de- 
scribing the graph we analyze 0-[l3 • The validity of all 
these approaches is still under investigation, see e.g. [ll|. 
Note however, that all these variations listed above work 
on equations describing relationships between probability 
vectors. In particular, for the last example, the system is 
being described by pi - probability, that i-th node is in- 
fected. The problem, however, is that there is not a single 
moment when a particular verte is - say - 0.41 infected. 



A vertex can be either infected (1) or not (0). This prob- 
lem has been already noticed by Petermann and De Los 
Rios [13. 

In this paper we introduce another model of epidemic 
spreading and analize it with complitely different ap- 
proach. Let us focus on sexual transmitted diseases. For 
this case the assumption that a particular node is able 
to contaminate more then one of its neighbots during a 
time step seems not to be the most suitable one. Bearing 
this idea in mind we develope a single infection epidemic 
spreading model. 

This paper is organized as follows. In section II we 
describe proposed model of epidemics spreading with at 
most one infection per time step. In Section III the theo- 
retical analysis of the model: epidemics treshold, station- 
ary state and mixing time, is being investigated. Simula- 
tions are presented in section IV. In Section V we draw 
the final considerations. 



II. MODEL DESCRIPTION 

Consider a connected, unweighted graph with n ver- 
tices enumerated by indices i = 1, . . . , n, described by 
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transition matrix {Pij}, ^.j^i Pij — 1- The model will 
be of SIS kind: all the individuals are at the beginning 
considered as susceptible (S). After contamination they 
become infected (I) but they still have a chance to recover 
and be susceptible again. 

We start thus with the all but one nodes susceptible. 
The one which is infected is chosen at random. The whole 
process consists of 4 actions which we repeat at each dis- 
crete time step. At each time step we choose randomly, 
with identical probability —, a node - say - i-th one. Then 
we choose its neighbor according to the transition matrix 
{Pij}, i.e. there is Pij chance that we point j-th vertex. 
If one of these two individuals J, j is infected, it contam- 
inates the second one with probability z. At the end of 
each time step we recover each infected node with prob- 
ability r. 

Note, that this method restricts not only each infected 
node to contaminate at most one of its neighbor. In fact 



we restrict all the dynamics to at most one contamination 
per time step. One can say it is not realistic approach. 
However, from the one hand we can say that this could 
be a basis for further generalisation. From the other: 
we can imagine, and even find in reality, systems that 
fulfill assumption described above. In large, academic 
cities there are often big flats situated in old tenement 
houses, settled by quite large amounts of students, who 
live with 3-4 roommates per chamber. As there is no 
space for privacy in this way of living, they sometimes 
devote one room in the flat to be a so-called sexroom, so 
contamination by sexually transmitted diseases can take 
place at most once per time step (say: per hour). This 
seems to be a good example of a system which can be 
described by our model. 



III. MODEL ANALYSIS 

In order to mathematically describe the model we de- 
fine Xj (t) which takes the value 1 if the node j is being 
contaminated by one of its neighbors at the time step 
t, and otherwise. Furthermore, we denote the set of 
all infected nodes at the time step t by I{t). We will 
be interested in the expectation value of Xj{t) with a 
condition that the set of infected nodes consists of some 
particular vertices. 

There are two independent ways of contaminating j- 
th node during one time step. Either we choose j-th 
node (with probability — ) and then one of its infected 
neighbor (with probability J2k£i(t) ^jk) or we choose j- 
th node's neighbor (with probability - for each one) and 
then we pick j'-th node (it happens with probability Pkj 
for a particular node k, so J2k£i(t) ^kj for all of them). 
Summing up we obtain: 



L{x,\m) = -[ J2 p^^ 

kei{t) 



kei(t) 



(1) 



where both terms are multiplied by contamination prob- 
ability z. We are, however, interested in the behaviour 
of whole system, not one node only. 

Let us thus define D{t) - expectation value of change 
of the number of infected nodes. Due to additivity of 
expectation value we can write: 

i^(i) = E(|/(m)i-i/(t)ii/(i))= J2 nx,\i{t))-r\m\, 

ant) 

(2) 

where, apart from adding all 'K{Xj\I{t)) terms, we sub- 
stract the term responsible for healing: number of in- 
fected nodes multiplied by recovery probability r. Using 
Eq. ([T]) we immediately conclude: 

^w = ^( E Pj^+ E Pk,)-r\m\. 

kei(t),jii(t) kei{t)j<fi{t) 

(3) 



The equation above defines the dynamics of the model: 
by solving it one could provide the complete informa- 
tion about the process. Unfortunately, in general sums 
T,keiit),jti{t)Pjk, T,k&i{t),j<^i(t)Pk3 are not precisely 
known as they strongly depend on the shape of the set 
I{t). We will show, however, that we are able to de- 
rive exact result for epidemic treshold for any graph and 
stationary state for some special cases. 



A. Epidemic treshold 

Our first aim is to find out the epidemic treshold for 
the process described above. We are interested in some 
relation of model parameters n, z, r that defines a bor- 
der between two situations: dropping and rising of the 
number of infected nodes in the beginning of the process. 

We are going to analyze Eq. ([3]). We have there two 
sums that look similar, so the first idea would be to 
add them somehow. But in general J2k£i(t) j4i{t) Pjk t^ 
X]feG7(t) iaiit)Pkj, so we cannot that easily simplify 
this equation (except for {Pij} - bistochastic). How- 
ever, in order to find epidemic treshold, we are inter- 
ested in the behaviour of the system in the vicinity 
oi t = 0. Let us thus take |/(i)| = 1 then, as it 
happens at the very beginning of the evolution, and 
denote the only infected neighbor by index I. Then 
J2kei{t),j<fi{t) Pkj = z2j(i{t) Pij ^ 1 (^'^ {Pij} - stochas- 
tic) and J2kei{t),3<fi{t)Pjk = Y.j0(t)P]i- Now we use 
the fact that at the beginning the first infected node is 
being chosen uniformly at random. Therefore the last 
term should be averaged over all possibilities of choos- 
ing ^: k T.iev T.j0(t) Pji = Ti T.j0(t) 1 = /'■«cn - In, 
where we used once again the fact, that {Pij} is stochas- 
tic. Finally we write the condition D{t) > which indi- 
cates the epidemic outbreak: 



D{t) < - X (1 



n-1 



) -r < 0. 



(4) 



For large n the term ^^i can safely be substituted by 
— . The condition for epidemic treshold for the model we 
proposed is thus: 



n 
2' 



(5) 



B. Stationary state 

Let us now turn to stationary state problem. The 
model being analized is by definition a purely Marko- 
vian one and above the epidemic treshold we anticipate 
our system to stay at some non-zero stationary state, i.e. 
we expect that the number of infected nodes will, in long 
times, oscillate about a fixed value. Practically however, 
due to statictical flucutation in finite real or simulational 
system, the epidemy may die out even above the treshold. 



The stationary fraction of infected nodes in general 
case (not specifying any particular shape of the graph) is 
not as easy reachable as the treshold calculated in the last 
section. What we basically have to do is to use once again 
all the formalism presented above and find the solution 
for the equation D{t) = without the constraint |/(i)| = 
1. The problem is to compute the sum X^fce/i't') nut) ^i^ 
- a task which is not trivial. We will thus estimate only 
stationary state for general case. In later subsections we 
give exact solutions for special cases of complete graph 
and uncorrelated homogenous graph. 

In order to perform estimation of the stationary state, 
we introduce the notion of graph conductance [l6| : 

Definition 1. Conductance of a given graph G described 
by a stochastic matrix {Pij} is: 



<I>(P) = min ■ 






scvmin{\Sl\V-S\}' 
where V is the set of vertices of a graph G. 



(6) 



C. Special cases 

1. Complete graph 

For complete graphs, i.e. graphs with all possible links 
present, we easily find the exact solution of stationary 
state problem. Note, that for this special case: 



\m\{n-\im 

n-l 



kei{t),j<^i{t) kei{t),j<^i(t) 

(11) 
as each of |/| infected nodes is linked to each of (rt — |/|) 
susceptible nodes by an edge chosen with probability ^^^y 
as each node has (n — 1) neighbors. We can thus find 
explicit and exact condition for D{t) = 0. From Eq. ^ 
we get: 



iit)s 



r{n — 1) 



(12) 



2. Uncorrelated homogenous graph 



This quantity measures how well-connected a given 
graph is. Due to the definition above we will analyze 
separately cases with the stationary fraction of infected 
nodes is = —f^ smaller and greater than i. 

Consider first is > \- Then also \Is\ > n — |/„| and, 
using Eq. ^, we lowerbound D(t): 



2z 
D{t) > — $(P)(n- |/|) -r\I\. 



(7) 



Bounding the latter expression in Ec^. ([7]) from zero we 
find that D{t) is positive for i < i < i, \n , therefore 



the stationary fraction is must be higher than this: 

1 



is > 



1 



(8) 



2z*(P) 



Let us now focus on the opposite case, mainly is < h, 
\Is\ < n~ |/„|. We again lowerbound D(t) using Eq. ([3|: 

D{t) > — $(P)|/| - r\I\ > — $(P)|/| - r(n - |/|).(9) 
n n 

Bounding right hand side of Eq. Q from zero, we con- 
clude analogically to the situation above: 



is < 



1 



1 + 



22$(P) 



(10) 



Let us consider now hypothetical uncorrelated homoge- 
nous graph. The term "uncorrelated" stands for the fea- 
ture that the probability that an edge departing from a 
vertex of degree j points on a vertex of degree k is inde- 
pendent from the degree of vertex j. By "homogenous" 
we mean that average number of connections between 
sets of vertices of some fixed sizes depends only on these 
sizes, not on the actual constituents of those sets. 

Bearing these assumptions in mind let us compute ex- 
pectation values of the two sums from Eq. ([3]) : 



= mECj:\k>l)\m\{n-\im, (13) 

where we put E(fc)/(n — 1) for the expectation value of 
existence of link between two vertices. We substract 1 
from n as a node cannot be connected with itself. The 
stationary infected nodes density comes to be: 



1 



r{n — 1) 



(14) 



where we denote (k) — K{k) and (1/fc) = E(l/fc). Specif- 
ically, for G{n,p) random graph (with the well-known 
binomial degree distribution) the product of {■^){k) goes 
to 1. In this case the latter result ([T4|) recovers the solu- 
tion for complete graphs ([T2l) . Moreover, ^(n,^) graphs 
are indeed uncorrelated in the limit of large n [13| , so we 
expect G{n,p) behaving like complete graphs for large n. 



This result, however mathematically correct, appears to 
be quite useless: the value of $(-P) is usually much lower 
than the sums that it approximates (X^fce/ft'i iai(t) ^jk, 
SfeG/(t) iii(t) ^kj) during the process. Let us thus work 
out exact results for some special cases. 



D. Mixing time 

In this chapter we will be interested in mixing time 
described in this article, i.e. the time needed by the pro- 
cess to reach the stationary state. Strictly speaking, this 



is kind of meta-stable stationary state, as in simulations 
on finite networks the only absorbing, stable state is the 
situation when the number of infected nodes is zero. It 
is clearly visible on Fig. [21 that we can distinguish two 
regimes with different behaviour: the regime of rapid in- 
crease in the number of infected nodes and the regime 
of stabilization. Let us state and prove a general the- 
orem restricting mixing time for any graph. The proof 
is inspired by related considerations for gossip spreading 
done by Shah fl^. 

Theorem 1. Let P be a stochastic transition matrix of 
a graph G of the size n. Then the mixing time T for the 
process described above fulfills: 

T{e) ^0{\ogn + \oge-^). 



Proof. We devide the proof into two parts, considering 
separately two stages of the process evolution: \I{t)\ < ^ 
and \I{t)\ > f . 

• \Iit)\ < f 

We recall first the general result for is < 5 stated in Eq. 
®: 

2z 
D(t) > — $(P)|/(t)| -r\I{t)\. 
n 

Denote now by A the smallest time t such that the num- 
ber of infected nodes exceeds ^ : 

A = inf{i:|/(t)|>^}, 



A A i = min(A,i). 

Note, that as long as \I{t)\ < f , we have A A (^ + 1) = 
A A t + 1. Recall now the general feature for any convex 
function g, Xi, X2 G R: 



l{xi) < g{x2) + g'(xi){xz - X2)- 



(15) 



Let us take: g{x) = -, xi = \I{t + 1)| and X2 — |^(i)|, 
then: 



< 



1 



^(|/(t+l)|-|/(t)|). (16) 



|7(i + l)| - |/(f)| \I{t + l)\ 
By construction of the process we have: 

\I{t + l)\<\I{t)\ + l^d\I{t)\, 

where 1 < d < 2, but as |/(^)| — 0{n) for n big enough 
the constant d can be arbitrarily close to 1. Now we 
continue with Eq. (J16p : 



where in the second line we used Eq. ^ and the defini- 
tion of D{t), Eq. ([2|). In the last line we used the fact 
that 1 — a; < exp(— z). Let us now define: 

cxp(ai) 



where a 



(18) 



Fk-3W)F(t*(^)I^WI-^I^Wl) ^ (17) 

^ (1 - (t$(P) - r)d--) < ^ exp(- J,(f $(P) - r)), 



We show that ({t) is a supermartingale, i.e. E(C(t)|{C(s) : 
s < t'}] < C{t) Vi' < t. As the only component of ({t) 
which is a random variable is I{t) and as the process we 
analyze is Markovian and asAA(t + l) = AAt + l, itis 
enough to show that E(C(A A (< + l))|/(A At)) < C(AAt). 
We do it using Eq. (fTTj) : 

E(C(AA(i + l))|/(AAt)) = 
exp((AAi)a)exp(a)E(^(^^Jj^^|/(AAt)) < (19) 

exp((A A t)a) exp(a) jjjxXtY^ exp(-a) = C(A A t). 

As ({t) is a supermartingale we conclude that E(C(A A 
t)) £ ]E(C(A A 0)) = 1. Furthermore, as we restrict our- 
selves to |/(i)| < f: 

2 
C(A At)> - exp((A A t)a), (20) 

n 

and directry from it we conclude that: 

n n 

E(exp((A A t)a)) < -E(C(A A t)) < -, (21) 

where in the last step we used the supermartingale prop- 
erty. Moreover, as exp((A Ai)a) f exp(Aa) as i — t- 00, we 
have also: 

n 

E(exp(Aa)) < -. (22) 

Finally, let us recall the Markov inequality: 

P(|X| > c) < ^^^ (23) 

c 

and choose ii = -{ln{n) — ln(e)). Then we straightfor- 
wardly get: 

™, . N ™/ /, N ^N E(exp(Aa)) e 

P(A > ti) = P(exp(Aa) > -) < ^ ); " < t^- (24) 

e 

• \m\ > f 

For this case we perform exactly the same procedure, 
but starting from Eq. ([7]) instead of Eq. (|9]), which we 
started with in the previous case. Following the same 
steps as above we only change constant a in Eq. (fT8)) 
into b = ^(|77<&(^) - ^^{P) - r), where explicitely 
appears the number of infected nodes at the stationary 
state. Second thing that has to be changed is Eq. ([20)) 
where, instead of ^ we can put n. Resulting time for this 
stage is: 



E(cxp(Ab)) , 

^ n _ 1 1 



P(A > t2) 

where ^2 = ■5(ln(?^) — ln(e)) 



(25) 
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Prom this general theorem we conclude, that the closer 
we are with chosen parameters to the zero-stationary 
state (i.e. the smaller is the stationary density of infected 
nodes), the slower is the first phase of rapid increase: 

Corollary 1. Mixing time is linear with inverse of the 
distance rj from the epidemics treshold, i.e.: 



r(77,e) = o(i(logn + loge-i) 



Proof. Recall Eq. ([3]): we demand D{t) > and trans- 
form this condition to: 



> 



\nt)\ 



kei(t).j<fi{t) kei{t)j<^i{t) 



(26) 



which boils down to equality for stationary state. We de- 
note right hand side of this equation by pc for the smallest 
possible situation, i.e. for epidemics treshold. Now let us 
take values of parameters z, n and r such that: 



Pc{l + v), 



(27) 



where r] > 0. Now we recall some parts of the proof 
of Theorem [TJ Actually, all we have to do is to rewrite 
condition for D(t) in parametrization given in Eq. (j27p 
and notion oi Pc- 



D{t)>r(--\m\~\I{t)\) ^rr,\m\. 
\nr Vr ' 



We put this result into Eq. pT|) obtaining: 



1 



1 t ^^ \ 

\iit+i)\-mT''^^^'' 



(28) 



(29) 



end then we proceed in the same way as in the proof of 
Theorem [TJ The result is 



^ E(cxp(Ag)) ^ ^ 



P(A > tc 



where tc = ^{Mn) - ln(e)). 



(30) 
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IV. SIMULATION 

Here we present simulations for stationary state of var- 
ious types of networks, i.e. complete graph, G{n,p) ran- 
dom graph [1^, Watts-Strogatz small world graph [17] 
and graphs with power law degree distribution (scale-free 
network, see e.g. [^)- Computer-simulational investiga- 
tions focus on the topics described theoretically in the 
last section, i.e. epidemic treshold, stationary state and 
mixing time. 




FIG. 1. Epidemic treshold for four different type of graphs: 
dots stay for simulational results, lines present theoretical 
prediction, Eq. ([5]). Starting from the bottom we have re- 
sults for G(n,p) random graph (green line, r = 0.006), scale- 
free network (bhie line, r = 0.008), complete graph (red line, 
r = 0.01) and small world graph (orange line, r = 0.012) 



A. Epidemic treshold 

We check here the behaviour of the process in the very 
beginning, i.e. exactly at the first time step. Four kinds 
of networks are being examined: complete graph, G{n,p) 
random graph with p = 0.5, small world graph with 
k — 6 neighbors on the circle and rewiring probability 
p = 0.5 (see [I4I) and scale- free network with the expo- 
nent a = 2.5. We vary sizes of networks n and for each 
type of the graph we choose different recovery probability 
r. Looking for the critical value of contamination prob- 
ability Zc we change parameter z and check for which 
value the fraction of infected nodes starts to increase. 
This procedure is being repeated 100000 times. Results 
are presented in Fig. [T] Visibly, simulations follow the 
theoretical prediction Eq. ([S]) prefectly for all four kinds 
of graphs being examined. 



B. Stationary state 

Results for stationary state are obtained by perform- 
ing many runs (typically 1000), finishing each of them 
at a fixed, long time step (10 000 - 100 000), cutting the 
beginning phase of rapid increase and fitting a line to the 
points oscillating about the stationary state. There are 
two types of results which we can end up with after a sin- 
gle run: epidemics either dies at a certain point (i.e. num- 
ber of infected nodes, due to fluctuations, reaches zero 
and - by construction of the model - stays zero, usually 
it happens at the very beginning of the process) or num- 
ber of infected nodes increases rapidely in the first stage, 
and then oscillates over some fixed value (see Figl^J. We 
call this value stationary state (presicely, as we have al- 
ready noted in Sec. IIIIDl meta-stable state). In order to 




FIG. 2. An example of a single run for z=l, r=0.005, random 
graph G(n,p) of the size 7i=100 and p=0.5. 




FIG. 4. Plot of stationary state value of infected nodes denstiy 
is for complete graphs versus contamination probability z: 
simulation (blue dots) and theoretical result (I12p (green line). 
We fix here n=100, r=0.005. 



0.5 -• 



FIG. 3. Plot of stationary state value of infected nodes denstiy 
is for complete graphs versus network size n: simulation (blue 
dots) and theoretical result (I12p (green line). We fix here z=l, 
r=0.005. 



compute average stationary state we neglect all the runs 
where there exist such a time step, when the number of 
infected nodes equals zero. 



FIG. 5. Plot of stationary state value of infected nodes den- 
stiy is for complete graph (blue circles), G{n,p) random graph 
with p — 0.1 (red squares), small world graph with k = 10 
neighbors on the circle and rewiring probability p — 0.5 (yel- 
low rotated squares) and scale-free network with the exponent 
a = 2.5 (green triangles) versus contamination probability z 
and fixed network size n = 100. 



First we examine complete graphs, as in the last section 
we provided the exact result for them p^ . In Fig. [3] we 
show how stationary infected nodes density is depends on 
network size n. Then, in Figj3J we show dependence on 
contamination probability z. Both figures show perfect 
agreement between simulation and theory, Eq. (J12p . 

As we have already seen the behaviour of complete 
graphs and how they relate to the theory described above, 
let us compare stationary state is for four different kinds 
of graphs. In Fig. [5] we show the results for complete 
graph, G{n,p) random graph with p = 0.1, small world 
graph with fc = 10 neighbors on the circle and rewiring 
probability p = 0.5 (see [iJl) and scale-free network with 
the exponent a = 2.5. Sizes of the graphs are fixed, n = 



100. Noticeably, the results for three out of four kinds 
of graphs are almost the same, while scale-free network 
goes an entirely different way. Below we will focus on 
complete, G{n,p) and small world graphs only. 



In Sec. nil C 21 we concluded, that G{n,p) graphs for 
large n should resemble like complete graphs. It is in- 
structive to see that in the limit of large n, epidemics, 
not only on G{n,p), but also on small world graphs be- 
haves the same as on complete graphs, see Fig. [S) 



C. Mixing time 

In this section we examine mixing times of the process, 
i.e. we check how long does it take to reach stationary 



...•^* 
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FIG. 6. Plot of stationary state value of infected nodes den- 
stiy is for G{n,p) random graph with p = 0.2 (blue circles) 
and small world graph with rewiring probability p = 0.5 (red 
squares) versus network size n. Number of neighbors on the 
circle k = 2n/10 is chosen such that the edges density ^ 
stays fixed. Red line shows theoretical prediction for com- 
plete graphs p2|l . We fix here z = \ and n x r = 1. 




FIG. 7. Average mixing time r for complete graph (blue 
dots), G{n,p) random graph withp — 0.2 (red squares), small 
world graph with rewiring probability p — 0.5 (yellow rotated 
squares) versus logarithm of network size ln{n). Number of 
neighbors on the circle k = 2n/10 is chosen such that the 
edges density ^ stays fixed. Simulational results are depicted 
by blue dots and red line shows theoretical prediction for com- 
plete graphs (fT2)) . We fix here r = 0.001 and n/z = 1000 in 
order to have stationary state not changed. Lines are plotted 
to guide the eye. 



state. Fig. [7] depicts how does average mixing time de- 
pend on ln{n), where n is network size, as usually. This 
is done for complete graph, G{n,p) random graph with 
p = 0.2 and small world graph with rewiring probability 
p = 0.5. For the same graphs we check average mixing 
time dependence on inverse of distance from epidemics 
treshold r] (see Corrolary [Ij . It is shown in Fig. [Sj 



FIG. 8. Average mixing time r for complete graph (blue 
dots), G(n,p) random graph withp = 0.2 (red squares), small 
world graph with rewiring probability p = 0.5 and fc = 20 
neighbors on the circle (yellow rotated squares) versus inverse 
of distance from epidemics treshold r). Simulational results are 
depicted by blue dots and red line shows theoretical prediction 
for complete graphs p2|) . We fix here r = 0.001 and n — 100. 
Lines are plotted to guide the eye. 



mixing time and show, that they are linear with ln(n) 
and l/ry, as theory in Sec. IIIIDI suggest by bounds of 
probability of mixing time proportional to ln{n) and l/i]. 



V. CONCLUSIONS 

We have proposed model of epidemics spreading with 
at most one infection per times step. Starting from the 
general formula for the change of the number of infected 
nodes ([3]) we provided condition for epidemics treshold 
for any kind of graph. Simulational results for epidemics 
treshold follow the theoretical predictions perfectly. Fur- 
ther more, stationary density of infected nodes for com- 
plete and uncorrelated homogenous graphs has been de- 
rived and bounds for this density, using the notion of 
graph conductance, have been obtained. Complete graph 
simulations show agreement with the theory. Epidemy 
on G{n,p) random graphs, according to no correlation in 
large n limit [l3|, as well as on small world graphs, in the 
large n limit, behave like epidemy on complete graphs. 

We have stated and proven theorem and corollary that 
bonds the probability of mixing time by values propor- 
tional to ln{n) and I/77, where n and 77 are size of the net- 
work and distance form epidemics treshold respectively. 
Simulations on complete, G{n,p) and small world graphs 
show even more, mainly that the average mixing time is 
linear with ln{n) and 1/ry. 
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